Interplay between structure and magnetism in M012S9I9 nanowires 
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We investigate the equilibrium geometry and electronic structure of M012S9I9 nanowires using ah 
tnitio Density Functional calculations. The skeleton of these unusually stable nanowires consists of 
rigid, functionalized Mo octahedra, connected by flexible, bi-stable sulphur bridges. This structural 
flexibility translates into a capability to stretch up to ~20% at almost no energy cost. The nanowires 
change from conductors to narrow-gap magnetic semiconductors in one of their structural isomers. 
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One of the major challenges in the emerging field of 
molecular electronics is to identify conducting nanostruc- 
tures with desired electronic properties, which are sta- 
ble and easy to manipulate. Due to their high stabil- 
ity and favorable electronic properties 0, carbon nan- 
otubes have been discussed extensively as promising 
candidates for molecular electronics applications. A ma- 
jor drawback is the current inability to synthesize or to 
isolate nanotubes with a given diameter and chiral angle, 
which determine their metallic or semiconducting nature. 
Moreover, single-wall carbon nanotubes are hard to iso- 
late from stable bundles, which form spontaneously dur- 
ing synthesis In this respect, recently synthesized 
nanowires based on Mo chalcogenides, such as M012S9I9, 
seem to offer several advantages by combining uniform 
metallic behavior, atomic-scale perfection and easy dis- 
persability 

Using ab initio density functional calculations, we 
study the suitability of Moi2SgIg nanowires as poten- 
tial building blocks of electronic nanocircuits. In terms 
of binding energy per atom, we find these nanowires to 
be almost as stable as carbon nanotubes. The nanowire 
skeleton consists of rigid, functionalized Mo octahedra, 
connected by flexible, bi-stable sulphur bridges. We find 
this structural flexibility to translate into an unusual ca- 
pability to stretch up to ~20% at virtually no energy cost, 
as a nanostructured counterpart of an accordion. Our 
calculations suggest that the nanowires change from con- 
ductors to narrow-gap magnetic semiconductors in one of 
their structural isomers. 

To gain insight into structural and electronic properties 
of these unusual systems, we optimized the geometry of 
infinite M012S9I9 nanowires as well as selected structural 
building blocks using density functional theory (DFT). 
We used the Perdew-Zunger form of the exchange- 
correlation functional in the local density approximation 
(LDA) to DFT, as implemented in the SIESTA code 
13 . The behavior of valence electrons was described b 
norm-conserving TrouUier-Martins pseudopotentials I, 
with partial core corrections in the Kleinman-Bylander 
factorized form We used a double-zeta basis, includ- 
ing initially unoccupied Mo5p orbitals. We arranged the 



nanowires on a tetragonal lattice, with one side of the 
unit cell given by the lattice constant a. The other sides 
of the unit cell were taken to be 19.7 A long, about twice 
the nanowire diameter, to limit inter-wire interaction. 
We sampled the rather short Brillouin zone of these ID 
structures by 6 k-points . The charge density and po- 
tentials were determined on a real-space grid with a mesh 
cutoff energy of 150 Ry, which was sufficient to achieve 
a total energy convergence of better than 2 meV/atom 
during the self-consistency iterations. 

Magnetic ordering in the nanostructures was inves- 
tigated using the Local Spin Density Approximation 
(LSDA) within the SIESTA code 7]. Even though 
LSDA may be considered less rigorous and dependable 
than LDA, we use it here as a systematic way to esti- 
mate net magnetic moments and exchange splitting in 
ID nanowires and finite segments of these systems. 
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FIG. 1: (Color online) Structural properties of M012S9I9 
nanowires. (a) Atomic arrangement within a unit cell. Planes 
normal to the wire axis are denoted by A, A', B, C, D. (b) 
Binding energy as a function of the lattice constant a. (c) 
Optimized inter-plane distances as a function of a. (d) Op- 
timum values of the Mo-S-Mo bond angle O as a function of 
a. Values in the two sulphur bridges per unit cell are distin- 
guished by different symbols. Dashed lines in (b-d) are guides 
to the eye. 
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With 30 atoms per unit cell, corresponding to 84 
degrees of freedom, global structure optimization is a 
formidable task. The initial experimental structure de- 
termination, based on Atomic-Resolution Transmission 
Electron Microscopy and X-ray diffraction data i pro- 
vided us with a useful starting point for determining the 
equilibrium structure. In spite of the significant energy 
gain of 50 eV with respect to the experimental structure 
during the initial structure optimization, the calculated 
atomic arrangement within the M012S9I9 unit cell, de- 
picted in Fig. n^a), was found to lie very close to the 
experimental structure. 

We found it useful to assign atoms to planes normal 
to the wire axis. We label these atomic layers A, A', B, 
C, and D in Fig.^a). Layers A and A' are structurally 
identical in terms of atomic arrangement, but rotated by 
180° with respect to each other. Mo octahedra, decorated 
by S and I atoms, form the structural motif of the AA' 
and CD bilayers. As we discuss later on, these bilayers 
can be thought of as rigid blocks, connected by bridges 
formed of three S atoms, denoted as layer B. 

Our results for the binding energy of the system as 
a function of the lattice constant a are summarized in 
Fig. n^b). For lattice constants other than the experi- 
mental value [ill aexpt — 11.97 A, the initial struc- 
tures used in the optimization were based on uniformly 
expanded or compressed nanowires, subject to random 
distortions, or structures optimized using semiempirical 
force fields. We considered a structure as optimized when 
different starting geometries resulted in the same struc- 
tural arrangement. 

The cohesive energy of the 30-atom unit cell, displayed 
in Fig.^b), translates into a formidable average binding 
energy exceeding 5 eV per atom. We found that inclu- 
sion of Mo5p orbitals provided a net stabilization of the 
system, but did not modify much the interatomic forces 
or equilibrium geometries. As indicated by our LSDA 
calculations, additional structure stabilization occurred 
due to magnetic ordering for selected geometries. 

Most intriguing in these results is the presence of mul- 
tiple structural minima. Since structures with optimum 
lattice constants a = 11.0 A, 12.3 A and 13.8 A are very 
close in energy, a distribution of bonding configurations 
along the chain axis will likely occur in order to maximize 
configurational entropy, resulting in the loss of long-range 
order. Identifying structures with very similar energies 
in a range of lattice constants also implies the possibility 
to stretch the nanowire by 20% at virtually no energy 
cost, similar to an accordion. The remarkable ability of 
MoeSj^Iy to expand easily by nearly 30% has been previ- 
ously noted in Atomic Force Microscope experiments on 
closely related systems 0|. This uncommon flexibility 
of an inorganic nanowire may explain in retrospect, why 
the observed lattice constant 0, 0| Uexpt — 11.97 A 
may differ from the structural minima in Fig. ^b). 

Trying to understand the origin of the unusual elastic 



behavior of Moi2S9lg nanowires, we first investigated the 
inter-layer spacings as a function of the lattice constant. 
Analysis of our data, presented in Fig.^c), suggests that 
interatomic spacing within AA' and CD bilayers changes 
very little with increasing lattice constant. The motif of 
the bilayers is formed of six Mo atoms in near-octahedral 
arrangement, with Mo-Mo bond lengths close to 2.6 A 
within the layers and 3.0 A along the wire axis. In the fol- 
lowing, we consider these decorated octahedra as rather 
rigid building blocks of the nanowires. We found the dec- 
oration order by I and S atoms to be non-random, since 
an interchange of these atoms within the unit cell raised 
the total energy by more than 0.5 eV. 

Most importantly, results in Fig. ^c) suggested that 
virtually all structural changes are accommodated by the 
sulfur bridges, which form the A'-D and C-A connections. 
We found that most changes occurred in the Mo-S-Mo 
bond angle 8, whereas the Mo-S bond lengths remained 
nearly constant. Our results for 8 as a function of the 
lattice constant a, shown in Fig.^d), suggest that the ex- 
pectation values of 8 do not change uniformly, but rather 
group around 90° and 150°. In the following, we will call 
the corresponding structural arrangements "short" and 
"long" bridges. Similar metastable structures have also 
been reported for ID wires of group IV elements In 
the case of isolated sulphur chains and S3 molecules, we 
have also observed a preference for similar bond angles 
as in Fig.^d). In analogy with similar preferential bond 
angles found in different allotropes of C and Si, we asso- 
ciate the short bridge with sp^ and the long bridge with 
sp^ hybridization. 

In the case of M012S9I9 nanowires, we find two such 
sulphur bridges in each unit cell. Local minima in the 
cohesive energy, plotted in Fig. Ub), are found at the 
lattice constants a = 11.0 A, 12.3 A, and 13.8 A. The fact 
that only three minima are found, and that a — 12.3 A 
is close to the average of the other two lattice constants, 
supports our conclusion that there are only two types 
of S-mediated bonds connecting the blocks, which are 
decoupled within the unit cell. The three local minima 
in Fig.^b) thus correspond to a short-short, short-long 
(or long-short), and a long- long bridge configuration. 

To verify that our optimum structures are not influ- 
enced by the constraints of a periodic lattice with a fixed 
unit cell size, we independently optimized the structure of 
the Mo-based building blocks with the proper decoration 
by S and I atoms, with sulphur bridge trimers attached 
at both sides. With the exception of the bridge atoms, 
which changed their positions, we found the atomic ar- 
rangement in MogSele, modelling the AA' bilayer, and 
M06S9I3, modelling the CD bilayer, to lie close to glob- 
ally optimized nanowire structures, discussed in Fig. ^ 

Comparing the binding energies of the bilayer clusters, 
E^ohiAA') = -81.0 eV and Ecoh{CD) = -89.2 eV, we 
conclude that the binding of sulphur atoms to the cluster 
is significantly stronger than that of iodine atoms. We 
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FIG. 2: Electronic band structures E{k) of M012S9I9 
nanowires with lattice constants (a) a = 11.0 A, (b) a = 
12.3 A, and (c) a = 13.8 A. Solid lines represent majority 
and dashed lines minority spin bands. Bands crossing the 
Fermi level are emphasized by heavy lines. Densities of states 
of nanowires with (d) a = 11.0 A, (e) a = 12.3 A, and (f) 
a = 13.8 A, convoluted with 0.01 eV wide Gaussians. The 
Fermi level lies at Ef = 0. 



also found that the isolated CD cluster acquires a net 
magnetic moment fi{CD) = 1.00 /ip. As we discuss in 
the following, not only finite nanostructures with dan- 
gling bonds, but also infinite nanowires may develop a 
net magnetic moment. 

We found the infinite M012S9I9 nanowires to gener- 
ally behave as conductors in the lattice constant range 
9 A< a < 15 A, with the exception of the metastable 
structure at a~ll A discussed below. In Fig. |21 we 
present results for the electronic structure at a = 11.0 A, 
12.3 A, and 13.8 A, corresponding to the equilibrium 
structures indicated in Fig. ^b). The interactions along 
the nanowire, causing band dispersion, depend crucially 
on the hybridization near the sulphur bridges. As dis- 
cussed earlier, changing the lattice constant does not af- 
fect interatomic distances, but rather modifies the Mo- 
S-Mo bond angle 0, thus changing the hybridization 
of directional orbitals. An increase in band dispersion, 
caused by increased hybridization, may be counter-acted 
by band repulsion in complex systems. 

The interplay between lattice constant and band struc- 
ture is depicted in Figs. [2Ia-c). Our results in Figs.E^b- 
c) suggest that the metallic character of the wires de- 
rives from a rather dispersionless band, which crosses the 
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FIG. 3: (Color online) Contour plots depicting the charge 
distribution in a Moi2S9lg nanowire at a = 11.0 A. (a) Total 
charge density contour at p = 0.05 el/ao''. (b) Difference 
charge density Ap(r), indicating regions of charge depletion 
and excess with respect to the superposition of isolated atoms, 
(c) p = l.OxlO"'^ el/ao^ charge density contour associated 
with states close to the Fermi level, Ef — 0.05 eV< E < 
Ef+0.05 eV. (d) Difference spin density p|(r)— p|(r), showing 
regions of excess majority spin t. Contour plots (b) and (d) 
are depicted in a plane, which contains the wire axis. 



Fermi level. As shown in Fig. I^^a), this band becomes 
very flat and even changes its slope at aw 11.0 A due 
to the changed hybridization associated with Mo-S-Mo 
bond bending. Consequently, the density of states de- 
velops a peak at Ep, which subsequently splits due to a 
magnetic instability, as shown in Fig.|3d). This splitting 
opens up a fundamental gap of 73 meV, accompanied by 
a magnetic moment of 1.00 /io per unit cell. At other 
lattice constants, the dispersion of this partly filled band 
increases, as seen in Figs. Elb) and (c). Consequently, 
as seen in Figs.Hfe) and (f), the density of states at the 
Fermi level is lower, suppressing the magnetic instability. 

A more detailed discussion of the electronic structure 
in the system with a net magnetic moment, at a = 11.0 A, 
is presented in Fig. |31 The charge distribution in the 
M012S9I9 nanowire, depicted in Fig.|3Ja), is very similar 
to the superposition of atomic charge densities. This is 
better seen in Fig. Olb), which displays the charge re- 
distribution within the system with respect to isolated 
atoms, given by Ap(r) = p(r) — Epatom(r)- As expected 
based on electronegativity differences, we observe a small 
amount of charge, stemming mostly from Mo4d orbitals, 
to be transferred to the region within the Mo octahe- 
dra, in particular to Mo-Mo bonds and sulphur atoms. 
The decorating iodine atoms do not experience a change 
in their charge distribution due to their presence in the 
nanowire structure. More interesting than information 
about all populated levels is the spatial distribution of 
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states close to the Fermi level, which form the valence and 
conduction band. Our results, shown in Fig. |2Ic), sug- 
gest that these states are rather delocalized. Together 
with the band structure results of Fig. |21 we conclude 
that the nanowires considered here should be conducting 
at room temperature, in agreement with experimental 
observations 

Spatially resolved information about the magnetic 
structure of the nanowire at a — 11.0 A is presented 
in Fig.|3[d). Associating up-spin with the majority and 
down-spin with the minority states, the plotted quantity 
p-|-(r) — px(r) allows us to identify spatial regions with a 
net magnetic moment. Based on the results in Fig. Efd) 
and a set of 3D spin density contours, not shown here, 
we find the majority spin to be localized mostly on the 
bridging sulphur atoms, in their non-bonding 2p orbitals 
normal to the wire axis, and to a smaller degree on ad- 
jacent Mo atoms. 

Since magnetism occurs in selected structural isomers 
of the infinite Moi2SgIg nanowires as well as in isolated 
building blocks, we expect the possibility of modifying 
the magnetic structure by locally stretching or expand- 
ing the nanowires. Since such structural changes may be 
induced locally at negligible energy cost, various mag- 
netic patterns could possibly be obtained by mechanical 
manipulation of M012S9I9 nanowires, possibly using an 
Atomic Force Microscope. Since magnetic ordering will 
likely affect spin transport in these systems, Moi2S9lg 
nanowires may be viewed as spin valves, which could 
change their state by applying local stress. 

We note that the overall metallic behavior of M012S9I9 
nanowires is rather unusual among chalcogenide struc- 
tures. Significant interest has been devoted in the past 
to the ability of layered semiconductors, such as MoSq, to 
form fuUerene- and nanotube-like structures |l,4 llfil Il7j . 
Similar to their bulk counterparts, M0S2 nanotubes were 
found to be semiconducting, unless intercalated with 
iodine to enhance conductance [11 [13. MoSz-based 
fullerenes and nanotubes, often containing iodine, were 
noted in particular for their favorable tribological prop- 
erties [isl ITtL i20j. At other stoichiometries, such as 
MogSale or M06S4I5, self-assembled nanowires have been 
reported to form with either semiconducting or semi- 
metallic behavior [21I E^ . In the particular M012S9I9 
stoichiometry studied here, the rich behavior of molybde- 
num chalcogenides is complemented by new properties, 
such as high structural flexibility, metal-semiconductor 
transition, and magnetism. We propose to use these sys- 
tems as unique building blocks in hierarchically assem- 
bled nanostructured electronic circuits. 
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